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Abstract In this paper, we explore the 3D structure of light scattering from dark-field illuminated 
live 3T3 cells in the presence of 40 nm gold nanomarkers. For this purpose, we use a high resolution 
holographic microscope combining the off-axis heterodyne geometry and the phase-shifting acquisition 
of the digital holograms. A comparative study of the 3D reconstructions of the scattered fields allows us 
to locate the gold markers which yield, contrarily to the cell structures, well defined bright scattering 
patterns that are not angularly tilted and clearly located along the optical axis (z). This characteri- 
zation is an unambiguous signature of the presence of gold biological nanomarkers, and validates the 
capability of digital holographic microscopy to discriminate them from background signals in live cells. 
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1 Introduction 



Gold nanoparticles attract great scientific and technological interest because of their physical and 
chemical characteristics. In particular, the optical tracking of gold nanoparticles in bi ology have gained 
popularity for several reasons. Gold nanoparticles provide high scattering efhciencies (jjain et al (20061)) 
and they can be detected directly using dark field or total internal r eflection (TIR) illumination down 
to particle diameters of 40 nm as shown bv iSonnichsen et~al ( 2000| ). Unli ke fluorescent ma rkers, they 



are immune to photo bleaching, and they are potentially non-cytotoxic f see I West et all ( 20061 )'). Because 
of these prop erties, the use of gold n anoparticl e s as b ioma rkers for live cell imaging using photothcrmal 
tracking (see ICognet et all (i20C)2') . Bover et al (|2003[ ) and lLasne et all ((200^)) has a high potential. 



As shown by Atlan et al ( 2008[ ). holographv has proved its ability to image and localize gold 



nanoparticles in 3D, either for flxed particles spin coated on a glass substrate or in free motion within 
a water suspension. More recently, lAbsil et a l (2010) have sho wn that heterodyne hologr aphy also 
allows the photothermal detection of 10 nm gold particles, and IWarnasooriva et"al (|2010[ ) have im- 
aged 40 nm gold particles in a cellular environment. In that last experiment the particle holographic 
signal is superimposed with the light scattered by the cell refractive index fluctuations, which yield 
a speckle fleld. For particles imaging, this induced speckle is a parasitic signal, but in many other 
situ ations, like in Dark Field micros copy , or in Differential Interference Contrast (DIG) microscopy 
(see iGoldberg and Burmeisteil ( 1986() ) this speckle is the main source of contrast that is used to image 



the cell itself. It is thus important to discriminate the particle signal from the cell parasitic speckle. 

In this paper, we have imaged biological s amples (3T3 cells) l abeled with 40 nm gold particles 
using the digital holographic setup described in IWarnasooriva et all ()201(]| ). We have performed the 3D 
holographic reconstruction of the wave-fleld scattered by the samples, which are illuminated at 45° 
in a total internal reflection conflguration, and we have shown that these wave-flelds are noticeably 
different for the particle, and the speckle signal. We showed here that important information can be 
derived not only from the intensity of the bright spots caused by the gold particles, but also from 
the 3D shape of the light scattering pattern, which is easily accessed using holography. We showed, in 
particular, that the speckle signal keeps memory of the illumination direction, while the particle signal 
does not. This result has been confirmed by imaging samples of cells that have not been labeled with 
gold particles, and samples of free gold particles. The shape of the wave-fleld scattered by the sample 
can thus be used as a signature that helps to discriminate the particle signal from the speckle. 



2 Materials and Methods 

2.1 Biological specimen preparation 

The biological specimens that we imaged are monolayers of live NIH 3T3 mouse fibroblasts labeled with 
40 nm gold particles via their integrin cellular surface receptors. Streptavidin-coated gold nanoparti- 
cles were attached to the surface cellular integrin receptors via biotin and fibronectin proteins: see 
FiglH Streptavidin and biotin are very well known for their strong affinity towards each other, and 
fibronectin, an extracellular matrix protein, has the property of interacting specifically with cellular 
surface receptors of integrin family. 

Fibronectinnroteins (fibronectin from bovine plasma, Sigma, St Louis, MO) were biotinylated 
using EZ-Link®Sulfo-NHS-LG-Biotin according to the provider protocol (Pierce, Rockford, IL). The 
final concentration of biotinylated-flbronectin solution was 0.447 mg/mL. The streptavidin-coated gold 
conjugates of 40 nm average diameter (BioAssay, Gentaur, France) were rinsed twice with IX PBS 
(Phosphate Buffered Saline) (pH = 7.25). We than diluted 10 of the gold solution in 990 /iL 
of the same PBS buffer solution. Then the dilute gold solution was incubated with 50 /iL of the 
biotinylated-flbronectin solution for four hours at room temperature to allow the specific streptavidin- 
biotin bonding. The final functionalized gold particles solution was kept at 4°G and used within 24 hours 
after its preparation in order to ensure maximum functionality. Before every use, the functionalized 
gold particles solution was sonicated. 
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Fig. 1 Gold bioconjugates and fibroblasts preparation procedures, (a) Gold bioconjugates functionalisation. 
(b) Fibroblasts-gold nanoparticles coupling. 

48 hours before the observation, monolayers of 3T3 cells were cultured in Duelbecco's modified 
Eagle's medium (DMEM Gibco, Invitrogen, Carlsbad, CA) supplemented with 10% fetal calf serum 
(PAA Laboratories GmbH) on 32 mm diameter fibronectin-coated glass cover slips (fibronectin from 
bovine plasma, Sigma, St Louis, MO) at 37°C in a 5% C02 atmosphere. After 24 hours of incubation, 
we added to each coverslip a solution composed of DMEM (2 mL) plus 500 fiL of the functionalized 
gold particles solution. The binding of integrin and fibronectin occurs at this level allowing the cells 
to attach, on their surface, the functionalised gold nanoparticles. 

The coverslip containing adherent 3T3 cells tagged with 40 nm gold nanoparticles was mounted 
on a specific observation chamber. In order to maintain the physiological pH condition during the 
experiments, cells were kept in DMEM-F12 medium (DMEM-F12 without Phenol red, B12 vitamin. 
Riboflavin, 0.5% fetal calf serum and supplemented with 20 niM of HEPES [(4-(2-hydroxyethyl)-l- 
piperazineethanesulfonicacid)] and L-Glutamine from PAA Laboratories). We measured the level of 
biotin incorporation on an HABA [2-(4' -Hydroxyazobenzene) Benzoic Acid)] quantitation assay to 
verify the efficiency of the biotinylation protocol. Average number of biotin molecules obtained per 
fibronectin is 2.5. 



2.2 Holographic Microscope Experimental Setup 

FigH] illustrates the optical setup. The illumination source is a single-mode near infrared laser diode 
emitting at A = 785nm (DL7140-201S 80 mW Laser Diode @90 mA current). A polarizing beam 
splitter cube (PBS) is used to split the original illumination laser light into two beams, a reference 
beam (complex field Eji, frequency fn) and an object illumination beam (complex field Eq, frequency 
/o) forming the two arms of a Mach-Zehnder interferometer. A combination of a half wave plate and 
two neutral density filters is used to prevent the saturation of the detector by controlling the optical 
power traveling in each arm. Two acousto-optic modulators (AOMl, A0M2) driven around 80 MHz 
and using the first order of diffraction, shift both frequencies at respectively fAOMi and fAOM2 ■ 

The object beam illuminates the sample, prepared as described in Section [2. 11 by provoking total 
internal reflection (TIR) at the medium-air interface in order to prevent direct illumination light from 
entering the system. The evanescent wave locally frustrated and the illumination wave directly scattered 
by the beads and cells give off a propagating scattered wave (complex flcld E) , which is collected by a 
microscopic objective (MO, 50 x magnification, NA=0.5, air). A beam splitter is then used to combine 
the scattered object wave and the reference wave which is slightly angularly tilted (~ 1°) with respect 
to the propagation axis of the object wave in an off-axis configuration. A half wave plate on the object 
illumination arm aligns the polarization of the corresponding beam ensuring its optimal interference 
with the reference beam. A CCD camera (Roper Cascade 512F, 512 x 512 square pixels of 16 fim 
size, exposure time 100 ms, frame rate fccD = 8 Hz) detects the interference pattern (hologram) and 
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Fig. 2 (a) Experimental setup : AOMl, A0M2: acousto-optic modulators; M: mirror; MO: microscope objec- 
tive (NA = 0.5); BS: beam splitter; PBS: polarizing beam splitter; CCD: CCD camera; Er'. reference field; Eo'- 
scattered field; z = 0: CCD plane; z = zq: CCD conjugate plane with respect to MO. (b) Details of the total 
internal reflection optical arrangement that is used for dark-field illumination of the biological sample. 

sends it to a computer. The hologram is then numerically treated and the complex field E{x,y, z) is 
reconstructed numerically. 



2.3 Holographic acquisition 

In order to filter out unwanted parasitic signals, we use a heterodyne modulation. A four-phase demod- 
ulation method is used to record holograms. This method consists in acquiring a sequence of images 
with a relative phase shift Atp — ■k/2 between two consecutive frames. The refere nce wave is frequenc y 
shifted by tuning the tw o acousto-optic mo dulators AOMl and A0M2 as done bv lLeClerc et all ( 200Clt ). 
and we get, as shown bv lAtlan et all (|2007^ ■ an accurate phase shift A(p that simplify the phase shifting 
digital holography demodulation. The heterodyne beat frequency is thus: 



^/ — fAOMl — fAOM2 



fcCD 



(1) 



where fccD =8Hz is the frame rate frequency of the CCD camera. The camera records a sequence of 
32 frames /o,...,/3i, and the object field E on the CCD plane (z = 0) is given by: 



i?(x,y,z = 0) = ^j "/„ 



(2) 



where ~ —1, and M = 32 is the number of frames used for the reconstruction. In Eqj^l the 
coordinates x,y (with < x,y < 511) are integers, which represent the pixel location within the CCD 
plane. The pixel size is then the physical CCD pixel size, i.e. 16 /xm. 



2.4 Holographic reconstruction 



The probl em of the reconstruc tion i n the context of holo graphic microsopy has been discussed in 
details bv IColomb et all ( 2006a|) and IColomb et all (l2006b). Neverthel ess, the Colomb et al. method 
refers implicitly to the phase-contrast imaging of iMargriet et a f (2005*), and is not well suited to the 
reconstruction of the 3D image of a wave-field as done bv lGrilli et al (,2001i) . Here, to get a 3D image of 
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the wave-field scattered by the gold particle, we have used a slightly d ifferent reconstruction method, 
which is inspired from the reconstruction method used by iMann eta 1 lloos). 

First, we considered that the measured hologram represents the field E{x,y,zo) within the CCD 
conjugate plane z = zq with respect to the Microscope Objective (MO), i.e., the plane, whose image 
is on focus on the CCD detector. Since we image the sample through a microscope objective MO, we 
must compensate the phase curvat ure, the phase tilt and the enlargement factor that are related to 
the presence of MO as disccused bv lColomb et all ( 2006bl) . We have thus: 



M 



E{x, y, z = zo) = e^i^-^+'^yy^ e^A(x^W) y j (3) 







where {Kx,Ky) and A are the tilt and lens parameters respectively that must be determined. 

We measured these parameters (or compensated their effect) by an original method that consists 
in reconstructing the image of the microscope objective output pupil by the one Fourier transform 
reconstruction method of Schnars and Jiip tncr (19M)- T he lens parameters A i s then close to the lens 
parameter that is used in the pupil reconstruction by the lSchnars and Jiiptned ( 1994) method. On the 



other hand the tilt parameters {Kx, Ky) are compensated by translating the pupil Schnars and Jiiptneil 
image in the center of the calculation grid. 
The properly compensated meas ured hologr a m rep resents then the field E{x,y,z — zq) in the 
conjugate plane zq. Then, as done bv lMann et all (|2005f ) in holographic microscopy, the field E{x, y, z) 
in the vicinity of the conjugate plane (i.e. for z ~ zq) is calculated by the angular spectrum meth od, 
which involves two Fourier transforms (see for example iLeClerc et all ( 200Clt ). ]LeClerc et all (|200l[ ) or 



lYu and Kim (2005)). This method is chosen here since it keeps the pixel size 5z constant whatever the 
reconstruction distance z is. 

The pixel size Sx = 6y, which must be calibrated to make a quantitative analysis of the holographic 
data, is measured by imaging with the same setup geometry a calibrated USAF target located in the 
CCD conjugate plane. We get Sx — Sy = 177 nm. The reconstruction is then done for 512 different 
reconstruction distances 

z = zo + (m^ - 256)(5z (4) 

where Sz — 177 nm and = 0...511. By this way, we get 3D volume images with 512 x 512 x 512 
voxels, with the same pixel size (6x ~ 5y = 5z ^ 177 nm) in the 3 directions x, y and z. 



3 Results and Discussion 



We have studied different samples of gold marked cell, unmarked cells and of free 40 nm gold particles. 



3.1 Samples of marked cells 



Figures 121 show the images of a first sample, with two cells and several particles. Figure |3Ua) shows a 
direct white light illumination image of the sample. The two cells can be seen, but the contrast is low. 
Figure [3] (b) shows the intensity holographic image of the same region of the sample reconstructed 
in the z = zq or z = 256 conjugate plane (here and in the following we will express the x, y and z 
coordinates by the corresponding pixel index mx^ my and m^). The display is made in logarithmic 
color scale. 

Because of the variation of the refractive index within the cells, the illumination light is diffracted 
yielding a speckle pattern that is superimposed with the particle signal. This speckle is visible on Fig. 
3 (b), and, from the envelope of the speckled zone, one can guess the shape of the cells. We interpret 
the brightest points 1, 2 and 3 of Fig. |2| (b), which correspond to maximum intensities / = 11.5 x 10^, 
4.6 X 10^, and 8.4 x 10^ and Digital Count (DC) respectively, as being particles signal. Many other 
bright points are also visible, but it is not simple to determine, which points are particles, and which 
are speckle hot spots. This is especially true within circle 4, where many bright points, close together, 
are visible. 

To better visualize the 40 nm gold pa rticles, we have displayed, on Fig.|3|(c), by using the Interactive 
3D Surface Plot plug- in of Image J fsee lAbramoff et al ((2004i )). a 3D linear surface plot of the region 
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Fig. 3 Images of the first sample with two fibroblast cells marked with 40 nm gold particles, (a) Direct white 
light illumination image, (b) Holographic intensity reconstructed image near the z = zO conjugate plane (z = 
255). The display is made in linear grey scale for the intensity I = \Ef. Black correspond to InJ = 5.9, white 
to In / = 16.3 Digital Counts (DC), (c) 3D linear surface plot of the same data. 



of the sample that is displayed on Fig|3](b). As seen, the particles that correspond to sharp peaks can 
be easily visualized, but some ambiguity remains around the meaning of the lower peaks, which could 
be attributed either to particles, noises, or scattering by biological features of the cell. 




Fig. 4 Images of the second sample with one fibroblast cells marked with one 40 nm gold particle, (a) Direct 
white light illumination image, (b) Holographic intensity reconstructed image for the z = 325. The display is 
made in logarithmic color scale for the intensity / = \Ef . Black correspond to In/ = 4.48, white to In/ = 15.56 
Digital Counts (DC), (c) 3D linear surface plot. 



The images of Fig|4]are obtained for a second sample with a single cell, and, as we will see, a single 
gold particle. Figure |3] (a) shows a white light image of the sample. One can see the cell, whose shape is 
triangular. FiguresH](b) shows the intensity holographic image reconstructed for z = 325 (i.e. 6.93 /im 
above the z = zq conjugate plane), with logarithmic color scale display. Here again, the speckle related 
to the light diffracted by the cell is visible on Fig|4] (b), and one can guess the triangular shape of the 
cell. The brightest point (arrow 1 on Fig. 2] (b)) is interpreted as a particle. Since the illumination 
intensity and focusing area is not well controlled, the particle maximum intensity / = 2.9 x 10^ DC 
obtained here is noticeably lower than for the first sample. Nevertheless, the signal obtained for the 
bright point marked by arrow 2 is more than 10 times lower (i.e. / = 1.9 x 10^), so we can interpret it 
as a speckle hot spot. This result is confirmed by Fig. H] (c) that shows a 3D linear surface plot of the 
sample. 

To go further, and to better characterize the particle's signal with respect to hot spots, we have 
analyzed the 3D images of the wave-field obtained by performing the holographic reconstruction for the 
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Fig. 5 Images of the second sample obtained by performing cuts of the 3D data parrallel to the y and z axis. 
Cuts are made in planes x = 128 (a) and x — 221 (2), which corresponds to the white dashed lines 1 and 2 of 
FigSJb) respectively. The display is made in the same logarithmic color scale than Figmb) by Volume Wiever. 



512 different reconstruction distances of Eq|31 The 3D data were displayed using the Volume Viewer 
plug-in of Image J, which is able to display cuts of 3D data in arbitrary planes. By using this plug-in, 
we have performed cuts parallel to the yz plane of incidence of the sample illumination beam (see Fig|2] 
(b)). 

Figure [5] (a) shows the image of a cut made in the plane x = 128, which intersects the particle 
signal 1 of Fig|3](b). The particle is seen as a bright spot on the yz cut image. We can notice here that 
the particle signal is located at 2 ~ 325 coordinate, which corresponds to the reconstruction plane of 
Fig|3](b). This is expected, since we have chosen to display on Fig|3](b) the plane where the maximum 
intensity is reach in 3D, and since this maximum correspond to the particle we consider here. 

The image of FiglS](b) corresponds to x = 221 i.e. to a cut that intersects the bright spot 2 of FigH] 
(b) we interpret as a speckle hot spot. Contrarily to the particle, the hot spot signal extension along 
the X axis (vertical axis on the Fig (5] images) is quite large. Moreover, the hot spot image is angularly 
tilted in the yz plane. 

This angular tilt can be simply interpreted by describing how the light propagation is governed in 
the biological cells. This propagation is dependant on the phase function inside the illuminated cell. 
Because biological tissues are inhomogeneous, the form of their phase function is not well defined and 
is thus characterized by the anisotropy coefficient g, which is the average cosine of the phase function. 
This parameter g describes the asymmetry of the single scattering pattern ; it is thus null when the 
scattering is isotropic, equals 1 for forward scattering and equals -1 in the case of backward scattering. 
In our experiments, the illuminated cells are maintained in DMEM medium that consists mainly of 
water. Since the refractive index of cells is close to that of water, the c ell anisotropy factor g is close 
to one {g ~ 0.9 in biological tissues as mentioned bv lCheong et all (119901 ) ). As a consequence, the light 
scattered by the cells mainly follows the forward scattering regime and the observed light scattering 
pattern appears to be tilted by approximately 45° since the incident illumination laser beam is initially 
tilted by 45° in conformity with the TIR illumination geometry (see Figl2](b)). 

The exact shape of the hot spot's wave-field can be calculated, but it is quite complicated. It involves 
the calculation of the angular distribution of the scattered light, which depends on the cell anisotropy 
factor g. One must then calculate the refraction of the scattered light on the medium-air interface, and 
take into consideration the collection of light by the microscope objective. The quantitative analysis 
of the wave-field's shape, which yields the angular tilt, is thus out of the scope in the present paper, 
and one can simply say that the hot spot signal keeps some memory of the illumination direction, and 
is thus angularly tilted in the yz plane. 

One can notice that a similar angular ti lt effect has been observed recently on the photothermal 
signal of 50 nm and 10 nm gold particles bv lAbsil et all ( 2010l) . 
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Fig. 6 Images of an unmarked cell sample, (a) Direct white light illumination image, (b) Holographic intensity 
reconstructed image for the z = 252. The display is made in logarithmic color scale for the intensity / = |-Ep. 
Black correspond to In J = 5.80, white to In/ = 13.28 Digital Counts (DC), (c) 3D linear surface plot, (d) 
Image obtained par performing cut of the 3D data in the plane x — 271 that corresponds to the white dashed 
line seen in (b). The display is made by Volume Viewer with the same logarithmic color scale as in (b). 

3.2 Control experiments performed on free particles or on unmarked-cells 

In order to confirm our interpretation of the angular tilt seen on FiglS] (b), we have performed some 
control experiments by imaging an unmarked cell sample and another sample of free gold particles. 

3.2.1 Unmarked cell sample 

FigE] shows the results obtained with a cell sample without particle. The direct white light image is 
shown in Fig. [BJa). Figure |6ljb) shows the holographic intensity reconstructed image that shows the 
light scattered by the cell refractive index inhomogeneities. This scattered light, which has a speckle 
structure, exhibits several hot spots that correspond to bright points on Fig[SJb), and we have chosen 
to image on FiglSJb) the reconstruction plane z — 252 that corresponds to the maximum hot spot 
intensity (5.85 x 10^ DC). 

The image we get here without particles is visually quite similar to Fig|3l^b) which is obtained with 
several particles. The 3D linear surface plot shown in Figl3jc) does not help since it is visually similar 
to FiglSJc). Nevertheless, the maximum of the bead signal (11.6 x 10® DC) in Figinijb) and (c) is much 
higher than the maximum of the hot spot signal (5.85 x 10^ DC) in FiglH^b) and (c). This makes the 
background noise visually higher in FigjS^c). 

Figure [S]Jd) shows the yz image obtained by performing a cut in the plane x — 271 that intersects 
the two brightest speckle hot spots of the sample. Now the yz image obtained in Fig. M,'^) without 
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Fig. 7 Images of free gold particles in a water and agarose suspension, (a) Holographic intensity reconstructed 
image for the z — 251. The display is made in logarithmic color scale for the intensity / — Black correspond 
to In / = 5.85, white to In / = 17.61 Digital Counts (DC), (b) 3D linear surface plot, (c) Image obtained par 
performing cut of the 3D data in the plane x = 258 that corresponds to the white dashed line seen in (a). The 
display is made by Volume Viewer with the same logarithmic color scale than for (a). 

particle is qualitatively different than with a particle in Fig. [SJa). The signal is angularly tilted like in 
Fig. Efb), and its extension in the z direction (vertical axis) is larger. 

3.2.2 Free particles sample 

The images of FigIT] are obtained for free particles (uncoupled to cells) in brownian motion in a water 
and agarose suspension. Here, we used agarose to slow down the motion of the particles enough to 
make the displacement of the particles negligible during the time of acquisition of the sequence of 
M = 32 frames. 

Figure IHfa) is the holographic intensity reconstructed image. The brightest points correspond to 
particles that are located within the reconstruction plane, while the less bright points are interpreted 
as out of focus particles. Reconstruction is made here in the plane z — 251 which corresponds to the 
maximum intensity (4.4 x 10^ DC) for the brightest spot (marked by a white arrow on FigEl^a)). Many 
particles can be seen on Fig. Hl^a) and on Fig. ^h) that displays the holographic data with 3D linear 
surface plot. 

Figure [SJc) shows the yz image in the plane x — 258 that intersects the brightest particle and 
corresponds to the white dashed line on Fig|6lja). As expected, the image of the particle (white arrow) 
has a smaller extension along the z direction (vertical direction) and does not exhibit the angular tilt 
we get with speckle hot spot in FiglSJ^b) and FiglHl^d). This result confirms our interpretation of the 
angular tilt observed with the speckle hot spots. 
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4 Conclusion 

To summarize, we have shown that heterodyne holographic microscopy, in the off-axis geometry, is 
well adapted to the detection of weakly scattering objects. The sensitivity, signal to noise ratio and 
selectivity of the technique allow the localization of light-scattering gold nanoparticlcs of a few tens of 
nanometers, which arc good candidates of non-toxic and perennial markers. Biological environments, 
however, are difficult to address since cell features generate strong parasitic speckle. Here, we have 
reported the detection of 40 nm particles attached to the surface of live 3T3 mouse fibroblasts. A 
comparison of these signals with either non-labelled cells or simple gold particles in solution allowed 
us to unambiguously discriminate particles. We show that, in addition to a stronger scattering signal, 
gold particles induce a relatively isotropic scattering, whereas biological features are characterized by 
mostly forward scattering. This dissimilarity in the scattering patterns, explained by the inconsistency 
of the refractive indexes and anisotropy parameters g, is easily characterized by digital holography, 
making it an excellent tool for the 3D detection of gold labels in biological environments. 

Acknowledgements Authors wish to acknowledge the French Agence Naionale de la Recherche (ANR) and 
the Centre de Competence NanoSciences lie de France (C'nano IdF) for their financial support. 
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